Plasmon mode modifies the elastic response of a nanoscale charge density wave system 
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The elastic response ol suspended NbSe3 nanowires is studied across the charge density wave 
phase transition. The nanoscale dimensions of the resonator lead to a large resonant frequency (~ 
10-100 MHz), bringing the excited phonon frequency in close proximity of the plasmon mode of the 
electronic condensate - a parameter window not accessible in bulk systems. The interaction between 
the phonon and plasmon modes strongly modifies the elastic properties at high frequencies. This is 
manifested in the nanomechanics of the system as a sharp peak in the temperature dependence of 
the elastic modulus (relative change of 12.8%) in the charge density wave phase. 



Research on nanoelectromechanical systems (NEMS) 
[TJ [2] has progressed a lot over the last decade, providing 
a way to explore the elastic properties of nanoscale sys- 
tems. However, the ability of NEMS to act as a probe 
of the lattice dynamics of a correlated system has not 
yet been widely recognized. This should be interest- 
ing because the screening of charge by electrons strongly 
affects the lattice stiffness in a metal. [3] The elastic 
properties are related to the dielectric constant - a fre- 
quency dependent quantity, and may be modified signif- 
icantly when the electrons undergo a phase transition to 
a coherent state. NEMS fashioned out of charge density 
wave (CDW) materials [2J [5] provide the right param- 
eter window to combine the advantages of a high me- 
chanical resonant frequency (in the MHz range) and an 
overlapping low frequency collective plasmon mode. The 
CDW systems studied in our experiment arc nanowires of 
the quasi-one-dimensional material NbSe3. Unlike bulk 
crystals, [6] the elastic modulus of these high frequency 
nanomechanical resonators shows an abrupt peak with 
the variation of temperature in the CDW phase. Using a 
model of an incommensurate CDW interacting with the 
ionic lattice, we attribute this surprising behaviour to a 
mutual locking of the lattice and the CDW. We discuss 
why this effect is visible only at the nanoscale. 

Crystals with a quasi-one-dimensional structure de- 
velop a charge density wave phase [H [5] at low tem- 
peratures where the electron density assumes a periodic 
variation over space along with a distortion of the lat- 
tice. The signatures of the new order parameter are ob- 
served in the electronic properties 18] and concomitant 
changes in the elastic properties [6l El [10] have been seen 
in macroscopic crystals. In recent years, electrical trans- 
port measurements on NbSe3 nanowires have shown that 
the properties of the CDW phase are significantly differ- 
ent from macroscopic samples. [TTHTS] The behaviour of 
elastic properties of nanoscale CDW systems are yet to be 
investigated in detail. Some measurements on thin flakes 
of the layered CDW material NbSe 2 have been reported. 





50 150 250 
Temperature (K) 




E 

h- 



65- 




DC Current (|iA) 

FIG. 1. (a), (b) Scanning electron microscope (SEM) images 
of suspended NbSe3 nanowires. (c) Four probe measurement 
of resistance as a function of temperature for a nanowire. 
(d) The CDW depins and participates in conduction under 
the application of a large electric field. This is observed by 
ramping up the dc current through the device and monitor- 
ing the differential resistance (dV/dl, where V is the voltage 
drop across the device and I is the current flowing through 
it). Above a certain critical current, the electric field is suffi- 
cient to depin the CDW and the differential resistance drops 
sharply. A lineplot at 46.8 K is shown (indicated by dotted 
line). The sharp drop indicates the depinning of the CDW. 
The threshold electric field for depinning is estimated to be 
11.4 V/cm. (Length of the nanowire = 3.1 /im.). 



P3] But the quasi-two-dimensional CDW phenomenon 
associated with this system is not well understood and 
a coherent picture of the underlying physics is missing. 
The CDW transition in the quasi-one-dimensional CDW 
system NbSe3 involves the gapping of a large part of the 
Fermi surface [TJ, leading to a very prominent influence 
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on its electrical transport properties. In this paper, we 
will probe the influence of the CDW on the elastic mod- 
ulus and mechanics of nanowires of NbSe3 . The coherent 
charge density wave has a considerable effective mass [TS] 
and combined with the large dielectric constant of NbSe 3 , 
[TO] reduces the CDW plasmon frequency low enough to 
leave its signature on nanomechanical measurements. 

Resonator devices were fabricated on a degenerately 
doped Si substrate (coated with a 300 nm dielectric 
layer of SiG^). The nanowires were suspended in dou- 
bly clamped geometry with metallic contacts at both 
ends [TT] (Figs. [l|a) and[ljb)). These were capacitively 
driven and set into mechanical vibrations to determine 
the resonant frequency. The change in elastic modulus 
was determined by measuring the resonant frequency of 
these devices across the charge density wave phase tran- 
sition. Typically, the width of the contacted nanowires 
are between 80-300 nm and the thickness is between 30- 
100 nm. The length of the suspended part between the 
contacts varies between 2-4 /im. The result of four-probe 
measurement of the resistance (R) as a function of tem- 
perature (T) curve of one such nanowire is shown in Fig. 
[ijc). The unit cell of NbSe 3 consists of three types of 
molecular chains, two of which undergo a CDW transi- 
tion. There are two CDW transitions observed in NbSe 3 
(marked T ci and T c2 in Fig. [TJc)). [7] 

Incommensurate CDWs (for which the CDW wave- 
length is not an integral multiple of the lattice spacing) 
[T8] can slide in the presence of an electric field. The 
CDW is usually pinned by impurities but beyond a cer- 
tain critical electric field (or a dc current passing through 
it), it depins and participates in electrical conduction. 
[8] IT9] Differential resistance measurements (Fig. [ljd)) 
on a nanowire device indicates the threshold electric field 
(E th ) for de-pinning of the CDW to be 11.4 V/cm at 46.8 
K. This is 3 orders of magnitude larger than macroscopic 
crystals. [8] A high depinning field as compared to bulk 
crystals has been reported to arise on the reduction of the 
size of the crystal. [T21 HO] The crossover lengthscale is 
determined by the phase coherence length [30] HI] of the 
CDW (~ few micrometers) . At low temperatures (below 
40 K), the resistance measurement plot in Fig. [ljc) shows 
metallic behaviour (resistance decreasing upon the reduc- 
tion of temperature). However, this is not the case for all 
our nanowires, and some of them show non-metallic [11] 
behaviour at low temperatures (see Fig. SI of Supple- 
mental Material [J7]). We now discuss the measurement 
of elastic properties by measuring the resonant frequency 
of suspended structures - the key aspect of our work. 

NEMS devices constitute an active field of modern day 
research and have been used in diverse applications like 
ultra-sensitive mass detection, [22], realizing bit-storage 
operations [23], and studying the quantum ground state 
of a harmonic oscillator [24] . The schematic diagram of 
our nanowire resonator device is shown in Fig. [2]ja). A 
frequency-modulated voltage Vf d M at frequency /, devi- 
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FIG. 2. (a) Schematic diagram of the circuit used for actu- 
ation and detection of mechanical resonance of a suspended 
nanowire. (b) Mixing current as a function of driving fre- 
quency and temperature for Device 1 (Vf c = -44 V). There 
is a large peak in resonant frequency (and also in the elastic 
modulus) below 60 K - where a charge density wave tran- 
sition takes place. (Inset: Lineplots of mixing current as a 
function of driving frequency at three different temperatures. 
The shift in the position of the peak (in 'frequency' axis) for 
the curve at temperature of 46.0 K compared to the other two 
curves corresponds to the peak seen in the main colourscale 
plot. The curves are plotted with an offset (along the 'mixing 
current' axis) for clarity.) (c) Variation of resonant frequency 
with temperature for Device 2 (Vg C = -45 V) . 



ation /a and modulation rate f r is applied at the source 
and a dc voltage Vg C is applied at the gate electrode. [25] 
(See Figs. S2 and S3 of Supplemental Material and the 
associated discussion about details of the measurement 
technique. [17) V s F d M (t) = Vbcos(27r/f + ^sin(27r/ r t)), 
where t denotes time and Vq is the amplitude of the ac 
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signal. This produces a driving force on the nanowire 
and sets it into oscillation at the same frequency /. 
The current through the device has a low frequency 
component at frequency f r (~ few hundred Hz) which 
shows a sharp change when the driving frequency / of 
Vf d M matches the natural resonant frequency of the 
nanowire. This component is called the mixing cur- 
rent I m ix and is measured using a lock-in amplifier. 

Imix oc ^/ (/2 ' / ° + ^ )(/2 - / "-^- ) where G is the con- 

ductance, q the capacitively induced charge, fo the nat- 
ural resonant frequency of the beam and Q is the quality 
factor of resonance. At the resonant frequency, the am- 
plitude of mechanical oscillation is maximum which leads 
to a sharp change of the mixing current. Therefore, by 
monitoring the mixing current as a function of frequency, 
we can determine the resonant frequency and from this 
the elastic modulus is inferred. The natural resonant fre- 
quency /o of a beam depends upon the elastic modulus 

E and strain r\ as / = \ff§ (gj \A?f + i s the 

length of the suspended nanowire, / the inertia moment, 
p the mass density and S is the cross-sectional area. [55]) 

To measure the change of elastic modulus at the CDW 
phase transition, the resonant frequency is determined 
by sweeping the frequency of Vf d M while measuring the 
mixing current, and the system is slowly heated up across 
the phase transition. The result from one device (Device 
1) is shown in the colourscale plot of I mix in Fig. [2^b) . 
The positions of sharp change in I m i X trace out the vari- 
ation of resonant frequency. There is a prominent peak 
(between 40-63 K) observed just below the T c2 CDW 
transition. Over the wide range of temperature from 
10 K to 190 K, the resonant frequency has a general 
trend of decreasing monotonically with increasing tem- 
perature, except that there is the abrupt peak between 
40-63K. The monotonically decreasing trend can be un- 
derstood by considering the strain in the system. When 
the device is fabricated at room temperature, strain in 
the system should be small. On cooling it down, the 
nanowire as well as the metal electrodes holding it con- 
tract. This tends to pull the nanowire and increases the 
strain and leads to higher resonant frequencies at low 
temperature, the change being a gradual one. By a sim- 
ilar argument, when the system is heated up starting at 
10 K, the nanowire and the electrodes will expand, thus 
reducing the strain and also the resonant frequency. This 
explains the gradual decrease of resonant frequency with 
increase in temperature. Now, in the region of the peak 
(at 40 K), the resonant frequency goes up abruptly. This 
can not be due to strain resulting from thermal expan- 
sion, because as explained earlier, strain would reduce 
with increasing temperature and will not lead to an in- 
crease in f . The quantity, a change in which can give rise 
to the 'peak', is the elastic modulus. Considering that 
/o oc VE, the height of the 'peak' corresponds to a frac- 



tional change in elastic modulus of 12.8% (±0.9%). This 
is an extremely large change. Comparison of the data 
of Fig. [2^b) with elastic modulus data on macroscopic 
crystals yields a striking contrast. In the vibrating reed 
experiments on bulk crystals by Brill et al., [B] there was 
no detectable change in elastic modulus at the T c2 CDW 
transition around 60 K. However, a 'dip' corresponding 
to a 0.09% change in elastic modulus was seen near the 
T cl transition at 142 K. We do not detect any significant 
feature around the T cl transition. (Our devices have 
Q~300 and this may not be sufficient for detecting the 
very small change in elastic modulus reported by Brill et 
al. [6] Further discussions are provided in Supplemental 
Material. [17] ) 

Fig. shows data from a second device (Device 

2). Two resonant modes can be noticed in this plot. 
Closely spaced modes arise in NEMS devices because of 
small asymmetries in device geometry. [27] For the upper 
mode, the signal of resonance becomes extremely faint as 
the resonant frequency starts increasing sharply beyond 
35 K and the 'peak' feature can not be traced (except for 
two points on the right hand 'slope' of the peak indicated 
by two arrows in the figure). However, the feature is clear 
for the lower mode and it is estimated that the maximum 
value of the 'peak' in resonant frequency corresponds to 
a 16.2 % change in the elastic modulus. (Data from one 
other device is presented in Fig. S4 of Supplemental Ma- 
terial. [T7] Measurements carried out on a total of six 
devices showed similar qualitative features.) 

In the following part of the paper, we discuss the mech- 
anism that leads to the abrupt variation of resonant fre- 
quency in our nanowire devices near the T C 2 transition, 
and also why such a feature was absent in measurements 
on bulk crystals. Phonons in metals arise from a Bohm- 
Staver mechanism where the restoring force for lattice 
displacements is provided by screening effects of the elec- 
tron gas. In a simple picture, an one-dimensional lattice 
can be thought of as a chain of masses connected by 
springs, each of stiffness K s and the stiffness scales in- 
versely with the dielectric constant e (K s cx -). [3] In 
the case of the CDW system NbSe3, the lattice also sees 
the potential due to the incommensurate charge density 
wave. The potential energy of the lattice, 

$ = \ K ^ x i ~ Zi-l " 2 + r (! ~ cos(2^)) . (1) 

i 

Here, T is proportional to the amplitude (or the or- 
der parameter) of the CDW, Xj is the position of the i th 
ion, I is the equilibrium separation between the ions and 
A is the wavelength of the CDW. (See Fig. [3] (a) for a 
sketch of the model.) This is known in literature as the 
Frenkel-Kontorova model. [28] In order to evaluate the 
resonant frequency, one needs to estimate the phonon 
modes. It had been shown by numerical computations 
[2"8] that when the ratio exceeds a critical value, a 



'phonon gap' opens up, shifting the acoustic phonon fre- 
quencies to larger values. Physically, the absence of low 
frequency modes at long wavelengths originates in the 
pinning of the lattice ions to the CDW potential. When 
the CDW potential is weak, the lattice vibrations are 
able to slide with respect to the CDW and one obtains 
the usual Bohm-Stavcr phonons. It had been observed 
in measurements of the dielectric constant of NbSe3 at 
microwave frequencies that e has a large peak around 
42 K (Fig. [37b), data reproduced from Gruner et al. 



[IB]). This implies that the parameter ^- (Fig. |3|c)) 
will have a peak in the similar temperature range - pre- 
sumably taking the system into the 'phonon gap' regime 
producing a large change in the experimentally measured 
resonant frequency. In Eq. 0, we have disregarded the 
vibrations of the CDW, assuming it to be of infinite mass. 
In reality, the CDW is coherent up to a length scale of a 
few microns and thus has a mass m* ~ 10 4 m e . |15j As 
this is of the order of the ionic masses that are associated 
with the phonons, we should also take into account the 
vibrations of the CDW as they will mix with the phonon 
modes. 

There are two ways in which the CDW can contribute 
to the apparent stiffness of the sound modes. The most 
natural one is that the opening up of the CDW gap in- 
creases the stiffness of the electron gas, which itself con- 
tributes to the elastic modulus assuming that the lattice 
and CDW are in equilibrium. In general, this is a small 
effect, of the order of (A/W) 2 ~ 1CT 6 , the squared ra- 
tio of the CDW gap (A) to the overall electronic band- 
width (W). This is consistent with measurements on 
bulk samples. A second route is more subtle, and is as- 
sociated with the fact that the CDW, when distorted by 
being pinned to the lattice, produces long range Coulomb 
forces, which are themselves screened by the metallic 
background in equilibrium. However, when driven at a 
frequency which is too rapid for the electron gas to re- 
spond, the frequency of a longitudinal sound mode will 
be shifted up to the 'CDW plasmon'. 

The sinusoidal charge distribution over space n(r) is 
represented as n(r) = n c + n cos[Q.r + 4>{r)). n is the 
amplitude and Q is the wave-vector of the CDW. n c de- 
notes the electron density in the metallic state. The spa- 
tial variation in the phase <p(r) of the CDW results from 
the deformation of the CDW due to pinning and applied 
electric fields. The dynamics of the charge density wave 
can be understood by analyzing the equation of motion 
of the CDW phase 4>(r). [15] The response of the system 
at frequency / is given by 
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FIG. 3. (a) The lattice is represented by a chain of masses 
connected by springs each with stiffness K a . The position 
of the i th ion is given by x%. The potential of the CDW 
is periodic with an amplitude T. K s ~ where e is the 
dielectric constant of the system, (b) Data of dielectric con- 
stant of NbSe3 at RF frequencies as a function of tempera- 
ture reported by Gruner et al. [16] (c) The parameter Y/K s 
is plotted as a function of temperature (T). It is evaluated 
following the relation K s ~ -, with the data in (b) for 45 
MHz frequency being used as a prototype for the variation of 
e with T. The potential due to the CDW is assumed to follow 
a Landau-Ginzburg expression: V ~ (T C 2 — T) 1 ^ 2 . 



o z and e z denote the conductivity and dielectric constant 
respectively along the chain direction (along which the 
CDW can slide). We have assumed for simplicity that 
this pinning is uniform (see Ref. [TSlfor a discussion). The 
last term of the above equation indicates that there exist 
two different frequency regimes where the dynamics of 
the CDW is distinctly different and a cross-over between 
these two regimes occur at a characteristic frequency of 
f c =<r z /(2TTe £z)- At low frequencies (/ < ^^), the 
response is overdamped, with an effective damping 7 = 
7o + n c e 2 /<j z . In this regime, the motion of the CDW 
is strongly screened by the conduction electrons. Any 
effect of the coupling of the CDW with the lattice (such 
as in the clastic properties) is hard to detect given the 
broadness of the response. At high frequencies, / 3> 
2nt t. ' t ne l as t term in Eq. (|2j) acts as a restoring term 
instead of a damping term, and one obtains a (CDW 
plasmon) resonance at a frequency f p given by 



Gl 1 = -4n 2 m*f 2 - i2n lQ f + V - i2nf ^— , 

(2) 

where, Vq is the CDW pinning potential coming from 
impurities and the lattice, 70 is the damping coefficient, 



f P = — v^o/m* + n c e 2 /m*e e z . (3) 

Using the conductivity data from electrical transport 
measurements on our devices and the dielectric constant 
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values as measured by Gruner et al., [T5] the crossover 
frequency f c = 2 ° e " Qe is estimated to be f c ~20 MHz. 
[17j This shows that our devices with resonant frequen- 
cies of tens of MHz exceed this crossover. In contrast, ex- 
periments for determining the elastic modulus of macro- 
scopic crystals by Brill et al. were performed with reso- 
nant frequencies not exceeding a few kHz (0.1-1.7 kHz). 
[6] It is not expected that the coupling of the CDW 
condensate to the mechanical motion of the ionic lattice 
would be discernible at such low frequencies. 

Now we estimate the CDW plasmon frequency f p . 
Assuming that the pinning potential can be neglected, 
f p ~ (l/27r)^n c e 2 /ra*e e z . We take n c ~2xl0 25 m" 3 
from an analysis presented in Ref. OHO The effective 
mass of the CDW m* -10 4 m e (Ref. QU). Using these 
numbers, we get f p ~39 MHz. 17 This is of the same 
order as the resonant frequencies of our nanowire devices, 
which means the Bohm-Staver phonons are close to the 
CDW plasmon frequency and the phonon modes can get 
renormalized due to a coupling between the ionic lattice 
and the CDW. It has to be kept in mind that the es- 
timate f p ~39 MHz is an approximate one. The value 
of n c used for the calculation is an order-of-magnitude 
estimate. [H] We have taken m e to be the bare elec- 
tron mass. It should be the effective mass of electrons in 
NbSe3, and was determined to be of the order of unity 
from magnetotransport measurements [SU]- The exact 
value of f p may be larger (though the order of magnitude 
should be the same) than the estimate provided. Thus, 
upon the renormalization of the phonon modes due to the 
closeness to the CDW plasmon mode, the elastic modu- 
lus of the NbSe3 nanowires can undergo a sharp increase 
as has been experimentally seen in our devices. 

On the basis of the Frenkel-Kontorova [35J model and 
the measurements of the high frequency dielectric con- 
stant by Gruner et al. [TB] , the following picture emerges 
for our experiments. For temperatures T » 42K, 
the resonance frequency corresponds to the longitudinal 
phonons corresponding to a wavelength of the order of 
the NbSe3 sample. In the Frenkel-Kontorova model, [25] 
this regime corresponds to free lattice vibrations unre- 
strained by the CDW (which is either weak or nonex- 
istent). At lower temperatures, when the CDW transi- 
tion takes place, the dielectric function is experimentally 
known to show a strong enhancement when measured 
at microwave frequencies of the order of 100MHz [IB] . 
This softens the Bohm-Staver phonons which now get 
locked to the CDW potential resulting in the opening of 
a phonon gap at long wavelengths. As the temperature 
is lowered further, the dielectric function again decreases 
[To] and the lattice vibrations are no longer locked to the 
CDW. 

This Letter demonstrates that the elastic properties 
of nanoscale resonators fabricated using a charge den- 
sity wave material are strongly affected by the plasmon 
mode of the electronic condensate. This provides us with 



a new window for interesting explorations of the coupling 
between electronic order and the mechanics of nanos- 
tructures. There exist a number of correlated solid state 
systems with intriguing phase transition behaviours in- 
volving interactions between the electrons and the lat- 
tice - such as materials with a structural Mott transition 
[3"T] and vortex phases in the peak effect regime of su- 
perconducting NbSe2- [32] Probing at the nanoscale can 
also highlight the behaviour of individual phase domains. 
This is a less explored avenue and leaves a great scope 
for future work on NEMS. 
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